Self- shielding effect of a single phase liquid xenon detector 
for direct dark matter search 
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Abstract 

Liquid xenon is a suitable material for a dark matter search. For future large scale experiments, single phase detectors 
are attractive due to their simple configuration and scalability. However, in order to reduce backgrounds, they need to 
fully rely on liquid xenon's self-shielding property. A prototype detector was developed at Kamioka Observatory to 
establish vertex and energy reconstruction methods and to demonstrate the self-shielding power against gamma rays 
from outside of the detector Sufficient self-shielding power for future experiments was obtained. 
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1. Introduction 
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Recent results from cosmic microwave background 
(CMB), large scale structure (LSS) and type la super- 
novae observations have yielded a standard model of 
cosmology: a flat universe consisting of more than 70% 
dark energy and about 23% dark matter with the remain- 
der ordinary (baryonic) matter [I, 2, 3]. The widely dis- 
cussed (nonbaryonic) dark matter candidate is Weakly 
Interacting Massive Particles (WIMPs). The newest and 
lightest particle, the neutralino (predicted by Supersym- 
metry (SUSY)), would be an ideal WIMP candidate fl. 

WIMP dark matter can elastically scatter off" the tar- 
get nuclei of a terrestrial detector, transferring energy to 
the nuclei. Dark matter can be observed by detecting 
this released interaction energy. But the main challenge 
is to identify a WIMP recoil event over many back- 
ground events. Recent best limits on an upper bound 
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for the cross section of WIMP spin independent interac- 
tions are given by the CDMSII |5], XENONIO [6] and 
ZEPLIN-III ITJ experiments. To increase the sensitiv- 
ity, a detector with larger mass and low background is 
required. 

XMASS was originally proposed as a 10-ton multi- 
purpose single phase liquid xenon detector aiming to 
detect pp-solar neutrinos, dark matter, and double beta 
decay [8]. For a large scale experiment, single phase de- 
tectors are attractive due to their simple detector config- 
uration and scalability. We have adopted a step by step 
approach and will initially develop an 800kg-sized de- 
tector (~ 100kg fiducial mass) dedicated to a dark matter 
search. Liquid xenon has a large light yield compara- 
ble to that of Nal(Tl) therefore allowing us to set a low 
energy threshold and achieve good acceptance for dark 
matter. Because of liquid xenon's large atomic number 
(Z - 54) and high density (~ 3 g/cnr'), background y- 
rays are quickly attenuated near the outer edge of the 
liquid xenon volume which makes a low background 
volume around the center of the detector. This self- 
shielding of the environmental gamma-Tay background 
is the key concept of the detector 

In this paper, we describe the results from a prototype 
detector by which we have measured basic properties 
of a single phase liquid xenon detector. We have es- 
tablished vertex and energy reconstruction methods and 
demonstrated the self-shielding power for background 
reduction. 



2. Experimental set-up 

2.1. Environmental backgrounds 

The WIMPs signal can be identified over many 
background events. For nuclear recoil events, the 
largest background contribution originates from neu- 
trons. Kamioka Underground Laboratory (the site for 
XMASS) is at a depth of 2700 m.w.e. (meters wa- 
ter equivalent), and has a smaller neutron flux than a 
surface laboratory by more than two orders of magni- 
tude The muon flux is also reduced by about five 
orders of magnitude |10] making the cosmogenic ac- 
tivation of xenon negligible. During the future opera- 
tion of the 800kg detector, environmental neutrons and 
gamma-rays will be reduced to negligible levels by a 
water shield of more than 2m thickness. However, back- 
grounds from PMTs and other materials used for the de- 
tector and supporting structures will not be reduced by 
such means mentioned above. Self-shielding is the most 
effective suppressor for those backgrounds. 



2.2. Liquid xenon scintillator 

The physical properties of liquid xenon are summa- 
rized in Table [T] 



Parameter 


Value 


Atomic number 


54 


Mass number 


13L29 [11] 


Boiling point at 1 atm 


165.1 Kfll] 


Melting point at 1 atm 


161.4K fll] 


Density at 16L5 K 


2.96 g/cm^' [12] 


Radiation length 


28.7 mm [13] 


Peak emission wavelength 


178 nm [14, 15] 


Spectral width (FWHM) 


~14nm[14, 15] 


Average energy required 
for one photon production 


23.7+2.4 eV [16] 
14.2 eV [17] 


Scint. Absorption length 


>100cm[18] 


Rayleigh scattering length 


60 cm [19] 


Refractive index 


1.61+0.1 [20] 



Table 1 : Physical properties of liquid xenon 

Existing reports vary on the average energy required 
to produce one photon by ionizing radiation. Also, the 
absorption length changes with the amount of water 
contamination in liquid xenon. We have checked these 
properties by using the prototype detector in order to 
demonstrate the feasibility of achieving the target sensi- 
tivity for a WIMP search with a large scale single phase 
detector. 

The relative scintillation efficiency of nuclear re- 
coils to electron recoils is called the quenching factor, 
/q. The nuclear recoil energy, £«, is converted into 
electron-equivalent visible energy, Ey, as follows: 

Ev = IqEr. (1) 

The unit of visible energy is keV electron equivalent 
(keVee), and the quenching factor of liquid xenon has 
been measured by several groups l,21i..22.,23i] . 

2.3. Detector 

The detector is a cubic vessel of 30cm in height 
made from oxygen free high conductivity (OFHC) cop- 
per with an inner volume of 27 litters. 

High purity 5N (> 99.999%) aluminum and 4N (> 
99.99%) MgF2 were evaporated on the inner surface as 
a reflective mirror for the scintillation light. 

Nine PMTs are attached to each exterior side of the 
vessel, and scintillation photons are detected by 54 
PMTs through MgF2 windows. Photo-coverage of the 



detector is about 16 %. Gaps between the PMTs and the 
MgF2 windows were filled with a high refractive index 
material, "Krytox® 16350," manufactured by DuPont (n 
= 1 .35 for 178 nm photons), allowing us to successfully 
increase the photoelectron yield of the detector by re- 
ducing total reflection between the MgF2 windows and 
vacuum. A schematic view and pictures of the prototype 
detector are shown in Figure s[T] and |2] 
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Figure 2: Pictures of the prototype detector. Reflectors are attached 
on the inner wall of the cubic vessel (top). Nine PMTs are attached on 
each side of the vessel and detect the scintillation photons (bottom). 



Figure 1 : Schematic view of the prototype detector. The cubic vessel 
containing liquid xenon is made from OFHC copper with an inner vol- 
ume of 27 litters. Scintillation photons are detected by PMTs through 
MgF2 windows. 



2.5. Radiation shields 

As the expected count rate caused by WIMPs is ex- 
tremely low, radiation shields are inevitable for dark 
matter search experiments. A set of shields as listed in 
Table[3]was adopted. 



2.4. Low background PMTs 

We have developed a low background 2-inch PMT, 
R8778ASSY, with Hamamatsu Photonics K.K. The 
PMT has reasonably high quantum efficiency, ~30 %, 
for 178 nm photons and can be operated at 173 K. Com- 
ponents of the PMT were examined individually and se- 
lected to reduce radioactive backgrounds. Table|2]shows 
the contamination of the radioactive impurities mea- 
sured by a high purity Ge (HPGe) detector at Kamioka. 



U (mBq/PMT) 


18 + 2 


Th (mBq/PMT) 


6.9 + 1.3 


K (mBq/PMT) 


140 + 20 


"'Co (mBq/PMT) 


5.5 + 0.9 



Table 2: Radioactive impurities in the R8778ASSY PMT measured by 
a HPGe detector The uranium and thorium chains were assumed 
to be in radioactive equilibrium. 



Outer shields 



material 


thickness 


polyethylene 


15 cm 


boric acid 


5 cm 


lead 


15 cm 


EVOH sheet 




OFHC copper 


5 cm 


Inner shields 


Inner OFHC copper shield 


6.7 cm 



Table 3: Material of the radiation shields. Listed from outside-in. 

The polyethylene shield and the boric acid powder 
shield act as a neutron moderator and absorber, respec- 
tively. The lead acts as the shield against environmen- 
tal y-rays. The Ethylene Vinyl Alcohol (EVOH) sheets 
with a reduced radon air purge system reduces radon 
gas from the inside of the shield. The OFHC copper 
shield acts as a shield against bremsstrahlung y-rays and 
X rays from the lead shield and is closest to the PMTs. 



2.6. Operation 

The detector and xenon gas line were evacuated and 
baked for two weeks at a temperature of about 373 K. 
The de-gas rate inside the detector and the gas Hne (total 
volume ~ 57 liters) after baking was on the order of 1 0"^ 
Pa/second at room temperature. Xenon gas was passed 
through a SAES getter (Model PS4-MT3-R-1) before 
filling the detector to remove impurities which shorten 
the attenuation length in liquid xenon. The radioactive 
impurity krypton cannot be removed with the getter so 
it was removed by a distillation tower before operation 
ill. 

A Gifford-McMahon (GM) refrigerator was used to 
liquefy the xenon. It also kept the detector at 173 K 
after filling with liquid xenon. For thermal insulation, 
the detector was suspended inside a vacuum chamber. 
The pressure inside the vacuum chamber varied from 
10"' to 10"^ Pa depending upon the temperature of the 
detector (room temperature to 173 K). 

2. 7. Electronics and data acquisition system 

The output signal from each PMT was fed into ADCs 
(V792, CAEN), a 250MHz FADC (STR7515, BAS- 
TIAN Technology), and leading edge discriminators 
(V814, CAEN). The gain of the 54 PMTs was set to be 
8.25x10^. The threshold level of the discriminator was 
0.25 photoelectrons, and the timing of each PMT was 
recorded by TDC modules (Tristan/KEK Online, TKO 
standard). The detector trigger required a coincidence 
of at least three hits within 95 ns. Trigger timing and 
veto timing were recorded by a TRG(TRiGger) module 
( VME standard) i26ll . The live time of the measurement 
was calculated from the TRG information. 

3. Detector Simulation 

A detector simulator program was developed by us- 
ing the framework of the GEANT3 package [27]. In 
the simulation code, tracks of particles, scintillation pro- 
cesses, propagation of scintillation photons, and the re- 
sponse of PMTs are simulated. 

In connection with the propagation of charged parti- 
cles, scintillation photons are generated according to a 
linear energyjransfer (LET) dependence of the scintil- 
lation yield 12811 given by the equation: 

dLldE = [A(dEldx)l(\ + BidE/dx))] + 770, (2) 

where dL/dE is normalized to 1 in the limiting case of 
dE/dx — > 00. A, B, and 770 are adjustable parameters, 
where 770 is the scintillation yield at zero electric field 
in the Umit of zero LET. In this paper, the values that 



manifest the best fit for liquid xenon ll28ll are used as t/q, 
A, and B. 



A^B^ 1.22. 



(3) 
(4) 



For the propagation of scintillation photons in liquid 
xenon, Rayleigh scattering and absorption are consid- 
ered in our simulation code. At the boundary of liquid 
xenon and other transparent materials, scintillation pho- 
tons are refracted or reflected according to the Fresnel 
equations. Light reflection and absorption in detector 
materials, such as the surface of the PMTs and the mir- 
ror on the inner wall of the detector, are also simulated. 
The light yield and absorption length of liquid xenon as 
well as the reflectance at the mirror on the inner wall 
were tuned to reproduce y-ray calibration data as dis- 
cussed in Sectio riSTI 

For the simulation of the PMT output of a single pho- 
toelectron, we employed an observed one-photoelectron 
distribution as a probability distribution function of the 
PMT output. For the simulation of multiple photoelec- 
tron signals, the PMT output was derived by summing 
up the PMT outputs for single photoelectrons with the 
required number of times. 

4. Vertex and energy reconstruction 

As discussed in Section!!] the self-shielding of liquid 
xenon for y-ray backgrounds is the key idea of the de- 
tector. Therefore the vertex reconstruction is a crucial 
issue for the detector. 

The vertex and energy of each event in the prototype 
detector is reconstructed with the light yield distribu- 
tion of the 54 PMTs. Light yield patterns of the PMTs 
from scintillation light at a single vertex are estimated 
by the detector simulation in a 30 cm cubic lattice with 
steps every 3 mm (1030301 points). The vertex is re- 
constructed with the following likelihood of a Poisson 
distribution by comparing the detected light yield distri- 
bution with the simulated ones: 



lo: 



exp ji" 



(5) 



Here, n is the observed number of photoelectrons and 

is 

/i, - Fiix, y, z) X num. of generated photons, (6) 

where F,(x, y, z) is the probability of observing one pho- 
toelectron at PMT No. / when a single photon is gener- 
ated at vertex {x,y,z). "num. of generated photons" is 



calculated as 



num. of generated photons = 



(7) 



where pe(i) is the light yield detected by PMT No. /. 

The vertex with maximum likelihood is selected as 
the reconstructed vertex. The reconstructed energy is 
then calculated as 



E [keV] 



"num. of generated photons" 

^ photon 



(8) 



Here, N photon is the number of generated photons per 1 
keV measured with 662 keV -y-rays from a '^^Cs source. 

5. Results 

y-rays were used to test various aspects of the detec- 
tor These y-rays are collimated through the radiation 
shield as shown in Figure [3] 
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Figure 3: Relative locations of the y-ray collimators (A, B, C) and the 
prototype detector. 



5.1. Measurement of xenon properties 

Characteristics of liquid xenon (e.g. absorption 
length and scintillation light yield) and the reflectance of 
the inside wall of the detector were measured by the col- 
limated y-rays from '^^Cs and *"Co. Here, the Rayleigh 
scattering length of liquid xenon for 178 nm photons 
was set to be 60 cm 1 19]. Light yield distributions were 
compared with those of the simulation by using the 



method, and the best fit parameter set with systematic 
errors were obtained as follows: 

• Absorption length = 66±10 cm 

• Reflectance of the mirror = 52+4-6 % 

• Average energy required for one photon production 
= 17.8+0.9 eV (obtained for 662 keV y-rays) 

These parameters are used for the following analysis. 

The average energy required for one photon produc- 
tion obtained for 662 keV y-rays is smaller than that es- 
timated for 1 MeV electrons by T. Doke et al. [ 161 and 
is larger than that obtained for an electron beam with 
kinetic energy < 100 keV by J. Seguinot et al. iIitIi . 

5.2. Event reconstruction 

The vertex and energy reconstruction methods were 
evaluated by using the y-ray data ('^'^Cs at 662 keV) 
from Collimator A, B and C. The X-Y projections of 
the reconstructed vertices of real data cluster around the 
y-ray injection points and agree well with that of the 
detector simulation as shown in Figure H) This demon- 
strates that the vertex reconstruction works as expected. 
The reconstructed energy spectrum of real data in suc- 
cessive fiducial volumes (full volume, 20 cm and 10 cu- 
bics from the center of the detector) agree well with that 
of the detector simulation as shown in Figure |5] and Ta- 
ble |4] Hence the energy reconstruction also works as 
expected. 
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Figure 4: The X-Y projection of the reconstructed vertices of Cs y- 
rays from Collimator A, B and C. The injection positions of the y-rays 
are shown in the bottom-right drawing. 
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Figure 5: Reconstructed energy spectrum of '^'Cs y-rays from Colli- 
mator A. Events reconstructed in full volume, 20 cm and 10 cm cubes 
from the center of the detector correspond to solid lines, dashed lines 
and broken lines. Photoelectric absorption peaks (662 keV) in each 
fiducial volume are fit with an asymmetric gaussian. The energy res- 
olution in each fiducial volume is shown in Table|4] 



5.3. Non-linearity of scintillation yield 

The reconstructed energy spectra of the following y- 
ray sources are compared with that of the simulation: 

• '^'"'Xe uniformly inside the detector (164 keV). 

• '^^Cs source from Collimator A (662 keV). 

• "^°K source inside the radiation shield (1461 keV). 

• ^°^Th source inside the radiation shield (2615 
keV). 

Events reconstructed in the 20 cm cubic volume are 
used for comparison. As shown in Figures |6] and |7] and 
Table |5] real data have a non-linear scintillation yield 
and energy resolution which cannot be explained by 
Poisson statistics. By introducing a LET dependence of 
the scintillation yield (discussed in Section|3]l, the simu- 
lation successfully reproduces the real data well except 
for ^°**Th. To reproduce the -'"^Th data with the sim- 
ulation, the scintillation yield must have a larger LET 
dependence than one expressed in Equation (|2]i. 

5.4. Self-shielding power of liquid xenon 

The self-shielding power of liquid xenon was evalu- 
ated with 7-ray data, '"Cs (662 keV) and ^°Co (1173 





mean (keV) 


resolution (%) 


Real data 


Full vol. 
20cm cubic 
10cm cubic 


670.9 + 0.7 
672.4+ 1.1 
676.1 + 1.8 


5.57 + 0.09 
5.47 + 0.13 
5.20 + 0.21 


Simulation 


Full vol. 
20cm cubic 
10cm cubic 


662.0 + 0.3 
663.4+ 1.1 
663.8 + 1.9 


5.83 + 0.07 
5.71+0.11 
5.59 + 0.20 



Table 4: Energy resolution of the prototype detector in each fiducial 
volume for 662 keV y-rays ('^'Cs from Collimator A). 
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Figure 6: Relative scintillation yield of real data (closed circle) and 
simulation (open circle) for various y rays. Vertical axis is normalized 
with 662-keV y rays ('^^Cs). Both real data and simulation have a 
non-linearity, ~ ±10 % from 662 keV. 



and 1333 keV) from Collimator A. If a y-ray vertex is 
near a PMT, it is difficult to accurately reconstruct the 
vertex because the ADC of the PMT saturates at 200 
photoelectrons. Therefore these ADC saturated events 
were removed. The reconstructed vertices along the y- 
ray incident direction (reconstructed z) was checked for 
events where the reconstructed vertices were within 10 
cm from the center on the xy plane as shown in Figure 
|8] Both ends of the distributions are distorted by cutting 
ADC saturated events. The number of events decreases 
as a function of the distance from the injected posi- 
tion of the y-rays which is expected from self-shielding. 
662 keV y-rays from '■'^Cs are attenuated by a factor 
of about 50 over a distance of 20 cm. While 1 173 and 
1333 keV y-rays from ^°Co are attenuated by a factor of 
about 3 over a distance of 10 cm. The distributions of 




Figure 7: Reconstructed energy resolution ot real data (closed circle) 
and simulation (open circle) for various y rays. 





mean (keV) 


cr(%) 


ijifflXe 


Real data 


183.5 + 2.0 


7.3 + 1.3 


164 keV 


Simulation 


180.3 + 0.2 


7.3+0.2 




Real data 


672.4+ 1.1 


5.5+0.1 


662 keV 


Simulation 


663.4+ 1.1 


5.7 + 0.1 




Real data 


1399 ± 15 


5.1 ± 1.1 


1461 keV 


Simulation 


1376 ± 12 


4.9 ± 0.7 




Real data 


2334 ± 15 


4.6 ± 0.6 


2615 keV 


Simulation 


2395 ± 18 


3.6 + 0.6 



Table 5: Fitted results of the photoelectric absorption peaks of various 
y-ray sources. Peaks are fitted with asymmetric gaussians. 
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Figure 8: Reconstructed vertex along the particle incident direction 
for 662 keV ('"Cs) and 1 173, 1333 keV (*Co) y rays. The distribu- 
tions of the real data (solid line) agree well with that of the detector 
simulation (dotted line). Both ends of the distributions are distorted 
by cutting ADC saturated events. 



real data agree well with that of the detector simulation. 



6. Conclusion 

A prototype single phase liquid xenon detector was 
developed. We have tested the performance of a single 
phase liquid xenon detector, have measured the physi- 
cal properties of liquid xenon, and have tested the ver- 
tex and energy reconstruction methods. By introducing 
a LET dependence in the scintillation yield, the simula- 
tion successfully reproduces the non-linear scintillation 
yield (+ 10 %) and the energy resolution of real data 
between 164 keV and 2615 keV. We have demonstrated 
that the self-shielding of liquid xenon for background 
y-rays is quite effective. 

The 800kg single phase liquid xenon detector is op- 
timized for a dark matter search. The photo coverage 



is ~ 70 % and the expected photoelectron yield is ~ 5 
p.e./keVee for an event at the detector center. This large 
photoelectron yield enables us to set the energy thresh- 
old to 5 keVee. With a fiducial volume cut using the 
reconstructed vertex, an ultra low background level as 
low as IQ^'^ I kg I day I keV can be achieved in the center 
40 cm diameter fiducial volume. The expected sensitiv- 
ity of the detector for dark matter is about 2 orders mag- 
nitude better than the present known limit |5] assuming 
5 years of data with 100 kg fiducial volume. 
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